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Is there a stellar upper mass limit? Recent statistical work seems to indicate 
that there is and that it is in the vicinity of 150 Mq. In this paper we use 
HST and ground-based data to investigate the brightest members of the cluster 
Pismis 24, one of which (Pismis 24-1) was previously inferred to have a mass 
greater than 200 M©, in apparent disagreement with that limit. We determine 
that Pismis 24-1 is composed of at least three objects, the resolved Pismis 24- 
ISW and the unresolved spectroscopic binary Pismis 24-lNE. The evolutionary 
zero-age masses of Pismis 24-lSW, the unresolved system Pismis 24-lNE, and 
the nearby star Pismis 24-17 are all ~ 100 M©, very large but under the stellar 
upper mass limit. 

Subject headings: binaries: close — stars: fundamental parameters — stars: 
individual (Pismis 24-lNW, Pismis 24-lSE, Pismis 24-17) — stars: early-type 



Introduction 



Massive stars strongly influence the structure and evolution of galaxies; thus, accurate 
knowledge of their fundamental parameters is critical. The spectroscopic binary frequency 
alone for OB stars is near 50 %; moreover, they are preferentially formed in compact multiple 
systems (IMason et al.lll998l ). Multiplicity therefore introduces a basic uncertainty in the 
determination of OB-star parameters, and it must be established with the greatest precision. 
Radial velocities are sensitive to separations up to a few AU; larger separations must be 
detected by means of high angular resolution. Even in th e near solar neighborhood, barely 



resolv ed or unresolved OB multiple systems are known (jWalborn et al.l 1 19991 : iNelan et al. 



20041 . and references therein). The unobservable gap in separations of course increases with 
distance, but it will be reduced with more powerful instruments. 



One important desiderat um is the stellar upper mass limit. The 03 (IWalbornlll97ll ) 
and 02 fjWalborn et al.ll2002l [W02]) spectral types include some of the most massive stars 
known. However, it is essential to know their multiplicities in order to determine their actual 
masses and solve the existing discrepan c y between the highes t accurate mass measured spec- 
troscopically (83 M0. lRauw et al.ll2004l : iBonanos et al.ll2004j ) and the photometric masses in 
the range 120-200 M© found by several authors. Another important desideratum is knowing 
how multiplicity influences the calculation of the initial mass function (IMF) at its high-mass 
end. Unresolved multiple systems artificially flatten the IMF by shifting objects from lower- 
mass bins into higher-mass ones. A sign of that effect is th e change in the log M — log L 
slope in the empirical relation derived from binary systems (INiemela fc GamenI 120041 ). At 
first glance, the Magellanic Clouds would appear to be the ideal place to study the fiattening 
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of the IMF, given their low extinction and foreground contamination and the small contribu- 
tion of depth effects. However, the very massive stars there are located typically an order of 
magnitude farther away than their Galactic counterparts; hence, for the Magellanic Clouds 
we are even more limited by the angular resolution of our instruments. 



Several authors flWeidner fc Kroupall2004l : lFigei]l2005l : lOev k Clarkell2005l : iKoenI 120061 ) 
have recently used the observed distribution of stellar masses in R136 and other young 
clusters and associations to establish that either (a) an stellar upper mass limit of ~150 M© 
exists or (b) the IMF becomes much steeper beyond 100 Mq. However, when W02 analyzed 
the earliest 0-type stars known, they found two apparent exceptions to that limit, the most 
egregious being Pismis 24-1 (also known as HDE 319718A and LSS 4142A). Two diff erent 
mass estimates fo r that star, based on the evolutionary calibration of IVacca et al.l (Il996l ) and 
on the models by ISchaller et al. (|l9£2|), yield 291 M^ a nd 210 Mq, re spectively. Pismis 24-1 
was c lassified as an 04 (f) star by lLortet et al.l (119841 ). as 03 HI by IVijapurkar fc Drilling 
Jl993h . as 03 If* bv lMassev et~alTj200lh . and as 03.5 If* by W02. Another star in the same 
cluster, Pismis 24-17 (also known as HDE 319718B), whi ch also appea r s to b e very massive 
fl Ol Ma, accordi n g to W 02), was classified as 03-4 V bv lLortet et aP fll984l ). as 03 III(f*) 
bvlMassev et al.l (120011) . and as 03.5 III(f*) by W02. Both objects were recently shown by 



Morrell et al 



(HVo have the same Oiv and Niv features in the 3400 A region as other 
stars with similar spectral classifications, thus confirming their spectral types. Pismis 24-1 
and 2 4-17 are the raain so urces of ionizing photons for the Hll region G353.2+0.9 in NGC 



6357 (IBohigas et al.l 120041 ). 



In this paper we present new observations that resolve the discrepancy between the 
~150 Mq stellar upper mass limit found in clusters and the larger apparent mass of Pismis 
24-1. 



2. Data 
2.1. Imaging 

In order to study the multiplicity at the high-mass end of the IMF, we were granted 
HST time in GO progr am 10602. The sample wa s drawn from an updated version of the 
Galactic O star catalog (IMaiz Apellaniz et al.ll2004l ) by choosing all stars with spectral types 
03.5 and earlier. A total of 20 objects with previous spectral classifications as 02-3.5 in 12 
fields were observed with ACS /WFC (the Wide Field Channel of the Advanced Camera for 
Surveys) in F435W+F550M-I-F850LP and with ACS/HRC (the High Resolution Channel of 
the same instrument) in a variety of filters (one of the objects had already been observed 
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in another of our programs, GO 10205). The Pismis 24 field was observed on 5 April 2006 
using the above mentioned filters for the WFC and F220W + F250W + F330W + F435W 
+ F550M + F658N + F850LP for the HRC (Fig. [I]). The HRC data were centered on the 
core of Pismis 24 and were dithered using 2- or 4-point patterns, depending on the filter. 
The individual HRC exposure times were selected to maximize S/N in all filters for the 
brightest stars while avoiding saturation. The resulting total integration times ranged from 
2 s (F850LP) to 704 s (F220W). The WFC data cover a large fraction of the cluster and 
were dithered using a 2-point pattern. The WFC exposures had total integration times of 
460 to 678 s for each of the three filters. No short WFC exposures could be added due to 
the large overheads involved; hence, the bright stars were heavily saturated, as evidenced by 
the obvious bleeding in the top panel of Fig. [H 

The ACS/HRC is an excellent detector for high-resolution imaging, with a pixel size 
of ~ 0'.'028 (the exact value depends on the position), a well-characterized PSF, an d a very 



precise geometric distortion solution (lAnderson fc Kingll2004j : iMaiz Apellanizll2007al ). Those 



properties allowed us to easily resolve Pismis 24-1 into two components of similar brightness 
(see bottom panel of Fig. [1]) in all of the seven HRC filterfl Given the similar magnitudes 
of the two components (see below), we will use the notation Pismis 24-lNE and Pismis 24- 
ISW in order to avoid possible confusions. Pismis 24-17, on the other hand, appears to be 
unresolved in all of the HRC filters, with its profile being essentially indistinguishable (apart 
from S /N effects) from that of the other bright isolated star in the HRC field. 

One of us (J.M.A.) has recently written a crowded-field photometry package, JMAPHOT, 
especially tailored for the HST instruments. For ACS, the code utilizes the SCLCOR (scien- 
tific, corrected) files which are the result of applying multidrizzle to the FLT (fiat-fielded) 
files in order to produce a DRZ (drizzled) frame. The SCLCOR files have the advantages of 
being free of cosmic rays and hot pixels and, at the same time, having the same geometric 
properties of the FLT files. Using geometrically uncorrected data allows for more precise 
position measurements and does not modify the photon statistics (i.e. adjacent pixels do 



not h ave correlated errors), thus improving the reliability of PSF fitting (lAnderson fc King 



20041 ). On the other hand, working on a geometrically distorted frame introduces some 
inconveniences, such as making the calculation of positions and aperture corrections less 
straightforward. When bright isolated stars are present, as is the case for our Pismis 24 
data, JMAPHOT can use them to model the PSF residuals and the aperture corrections. 



^The WFC, with a pixel size approximately double that of the HRC, can also resolve the system. The 
strong saturation prevents this from being immediately obvious in the top panel of Fig. [1] However, a closer 
inspection of the bleeding trail that originates from Pismis 24-1 reveals that it is double. 
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Fig. 1.— HST/ACS imaging of Pismis 24. (top) 2' x 2' FOV from the WFC F850LP drizzled 
data, (middle) 20" x 20" FOV from the HRC F850LP drizzled data, (bottom) 1" x 1" FOV 
from the HRC F435W drizzled data. In all cases N is towards the top and E towards the 
left. In the top and middle panels the intensity scale is linear while in the bottom one a 
logarithmic scale was selected to highlight the PSF. The bright stars in the top panel are 
heavily saturated. 
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We have applied JMAPHOT to the Pismis 24 HRC data. The results for Pismis 24- 
INE, 24-lSW, and 24-17 are shown in Table [TJ The magnitudes have been corrected for 
CTE (Charge Transfer Efficiency) and aperture effects and their uncertainties are calculated 
from the measured S/N plus a contribution from those effects. The aperture correction for 
F850LP requires a sp ecial treatment due to its strong dependence with the spectral energy 



distribution or SED (jSirianni et al.l 120051 ): we have used a value in accordance with the 
large extinction of Pismis 24 and we have also added 0.05 magnitudes in quadrature to the 
F850LP uncertainties to account for this effect. Also, as of the time of this writing there is 
no low-order flat field calculated for the F220W and F250W filters of ACS/HRC. For that 
reason, we have also added an additional factor of 0.05 magnitudes to the uncertainties for 
those filters. 

Our HRC data consist of 24 individual exposures. Combining the individual JMAPHOT 
measurements for the separation and position angle with their corresponding uncertainties 
we measured values of 363.86±0.22 mas and 207.816±0.015 degrees E of N, respectively. 



2.2. Spectroscopy 

2.2.1. Las Campanas 

Spectroscopic observations of Pismis 24-1 NE and SW were obtained with the Inamori 
Magellan Areal Camera and Spectrograph (IMACS) on the Magellan I (Baade) telescope 
at Las Campanas Observatory (LCO), on May 9.41 UT, under excellent seeing conditions 
(0.4" - 0.5"), which are characteristic of the Magellan telescopes. We used the f/4 (long) 
camera with a projected slit width of 2.08 A to cover the wavelength range 3660-6770 A at a 
reciprocal dispersion of 0.38 Apx^^. A spectrum of Pismis 24-17 was obtained a few minutes 
later, with the same instrumentation but somewhat poorer seeing (~ 0.6"). A He-Ne-Ar 
comparison spectrum and a fiux standard were also observed, as well as the usual series of 



Table 1. ACS/HRC photometry expressed in VEGAMAG 



Object 


F220W 


F250W 


F330W 


F435W 


F550M 


F658N 


F850LP 


Pismis 24-lNE 
Pismis 24-lSW 
Pismis 24-17 


16.628±0.051 
16.741±0.051 
17.504±0.051 


14.513±0.051 
14.621±0.051 
15.451±0.051 


13.103±0.010 
13.202±0.010 
14.027±0.012 


12.781±0.006 
12.897±0.006 
13.654±0.006 


10.996±0.017 
11.111±0.017 
11.783±0.023 


9.846±0.005 
10.030±0.005 
10.640±0.004 


8.360±0.050 
8.510±0.050 
9.011±0.050 
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bias and dome flats. These data were processed at LCO, making use of standard IRAI^ 
routines, to yield two dimensional wavelength and flux calibrated images. 

Given the separation of only slightly more than 3 pixels between the two Pismis 24 



1 com ponents in the IMACS Pismis 24-1 spectra, we used MULTISPEC ( iMafz Apellaniz 



20051 ) for the extraction. That IDL code is designed to obtain spectra from slitless or long- 
slit exposures of crowded fields using PSF fitting along the direction perpendicular to the 
dispersion. The extracted spectra for the three massive stars are shown in Fig. [21 



2.2.2. CASLEO 



Since iLortet et al.l (119841 ) had suggested from observed radial velocity variations that 
Pismis 24-1 probably is a spectroscopic binary, we acquired high resolution spectra to test this 
hypothesis. Spectral images of Pis24-1 were obtained with the 2.15-m J. Sahade telescope at 
Complejo Astronomico El Leoncito (CASLEO) in San Juan, Argentina, during two observing 
runs: September 1998 and June 2000. We used the REOSC echelle spectrograph with a 
TEK 1024x1024 pixels detector. The spectral resolving power was R ~ 15000. We obtained 
5 exposures covering the wavelength region 3800-6100 A and 3 exposures covering 4400- 
6700 A. We selected a slit width of 2" for all of our REOSC spectra, as this was the average 
seeing during our observations. Exposure times for the stellar images ranged between 30 and 
45 minutes, resulting in spectra of signal-to- noise ratio ~50. 

All the spectral images obtained at REOSC were processed and analyzed with standard 
IRAF routines. We determined radial velocities of the four most conspicuous absorption 
lines in our spectra, namely the hydrogen absorptions H7 and H/3, and the He ll absorptions 
4541 and 5411 A, fitting gaussian profiles within the IRAF routine splot. Also, nebular and 
interstellar lines were measured, and their velocities remained constant within the errors 
('2-3 km s'^) 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 2. — Rectified linear-intensity spectrograms of Pismis 24-lNE, 24-lSW, and 24-17 ob- 
tained witli tlie IMACS-Long Camera of the Baade telescope at Las Campanas Observatory. 
Each spectrogram has two segments, a violet one and a blue-green one, joined at 4370 A. 
The violet segments have been smoothed with 3-, 3-, and 5-pixel boxes, respectively (from 
top to bottom). The blue-green segment for Pismis 24-17 has been smoothed with a 3-pixel 
box. Constant values of 0.5 and 1.0 have been added to the Pismis 24-lSW and 24-lNE 
spectrograms, respectively. The non-dib lines identified are, in order of increasing wave- 
length, Hei+ii A4026, Niv A4058, Siiv AA4089-4116, Heii A4200, Hei A4471, Heii A4541, 
Nv AA4604-4620, Niii AA4634-4640-4642, Si IV A4654, Civ A4658, and Heii A4686, with dib 
standing for diffuse interstellar band. 
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3. The most massive stars in the core of Pismis 24 

3.1. Results from spectroscopy 

The spectral type of Pismis 24-lNE derived from the IMACS data, 03.5 If*, is consistent 
with that assigned by W02 to the composite observation of the system by Massey et al. 
(2001). Indeed, that spectrum defined the new type introduced by W02. The key criterion 
is the comparable intensity of N iv A4058 and N ill AA4634-4640-4642 (which stretches the 
original meaning of "f*" that the Niv is stronger than the Nlll). Pismis 24-lSW is of later 
type from both the N iv/N ill ratio (the former is weak but still detected) and the weak He I 
A4471 absorption line. Indeed, it now appears that the latter feature, as well as the strength 
of Hel+ll A4026, in the composite data is due to the SW component. Such is consistent 
with the nearly identical magnitudes of the two components and the relative morphology of 
their spectra. The spectral type of Pismis 24-lSW derived from the IMACS observation is 
04 III(f+), where the luminosity class follows from the weakened absorption and incipient 
emission at Hell A4686, and "f+" means that Si IV AA4089-4116 are in emission in addition 
to the Nlll (which is subsumed in the definition of f*). The spectral type of Pismis 24-17, 
03.5 III(f*), is also unchanged from the discussion by W02. The overlap between giant 
and supergiant absolute visual magnitudes at the earliest spectral types is discussed by 
W02; in Pismis 24-1 at least one additional component contributes to the relative apparent 
magnitudes of the resolved pair. 

Our REOSC spectra were obtained with worse seeing than the IMACS ones and it is 
obvious that they contain the combined spectrum of both visual components, Pismis 24- 
INE and Pismis 24-lSW, which are less than 0'.'4 apart. The absorption lines in our spectra 
appear to have variable shape, but we did not observe clearly separated component lines 
in our spectra. We determined the radial velocities of the mean profiles of the absorption 
lines, with the aim of verifyng the previously published radial velocity variations, which most 
probably would imply that at least one of the visual components of Pis24-1 also is a short 
period spectroscopic binary. The mean radial velocities of the absorption lines measured in 
our eight spectra present variations from -|-20 to —90 km s~^ (heliocentric) , thus confirming 



with somewhat larger amplitude the variations published by lLortet et al.l (119841 ). Since we 
are measuring the lines blended with the visual companion, the actual amplitude of the 
radial velocity variations may be considerably larger. Our data are too sparse to deduce a 
value for the orbital period, suffice it to say that it seems to be of only a few days. 

Our spectra show the signatures of a turbulent and high temperature interstellar medium 
commonly observed in H II regions surrounding very hot stars: namely that the interstellar 
lines have multiple components, and the presence of nebular absorption of He I 3888A. The 
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main components of the interstellar absorptions of Nal Dl and D2 in our spectra have 
heliocentric radial velocities of —13 ± 2 km s~^ and —12 ± 3 km s~^. 

Phil Massey and collaborators (private communication) have detected optical variability 
in the unresolved Pismis 24-lNE+SW photometry with a peak-to-peak amplitude of 0.07 
magnitudes and a period of 2.36088 days. This variability, together with the observed radial 
velocity variations, is a clear sign that Pismis 24-1 is at least a triple system. Furthermore, 
the radial velocity of the N iv A4058 line in the IMACS spectrum coincides with the value 
derived from our radial velocities obtained from the REOSC spectra at the appropriate 
orbital phaself]. Since N iv A4058 is visible only in the spectrum of Pismis 24-lNE, it must be 
the spectroscopic binary. We also point ou t that the fit of the X-ray spectrum of Pismis 24-1 



requires a two-temperature plasma model (jWang et al.ll2006l ). which is usually an indication 



of wind-wind interaction in a relatively close binary. The results here point towards Pismis 
24-lNE as the source of the hard X-ray component. 



3.2. Results from photometry 

The results in Tabled] show that Pismis 24-lNE is brighter than 24-lSW in the opti- 
cal/NUV by ~ 0.1 magnitudes, with nearly identical colors. Given their very similar spectral 
type and location, this is as expected: both stars have almost the same intrinsic colors and 
their light is affected by the same amount and type of dust. The most significant difference 
is in F550M-F658N [^-(Ha+continuum)], where Pismis 24-lNE is redder by 0.069±0.025 
magnitudes. This is also expected, since 24-lNE is a supergiant and has a stronger Ha 
emission (this effect is also seen in our LCO spectra). Pismis 24-lNE is brighter than 24-17 
by 0.924±0.016 magnitudes in F330W but the latter is slightly redder than either of the 24-1 
components; hence, the difference at F850LP is reduced to 0.651±0.007 magnitude^. The 
likely source for the color differences is that the light from Pismis 24-17 is more extincted 
than that from 24-1. 

We have processe d the Pismis 24-lNE, 24-lSW, and 24-17 photometry with CHORIZOS 



( iMafz Apellanial2004j ). a Bayesian code that compares observed magnitudes with SED fam- 



ilies and calculates the likelihood distribution. The current version of the code includes the 



most recent calibration for optical and NIR photometric zero points (iMaiz Apellanizll2007bl ). 



•^Note, however, that we have insufficient points to determine a fuU velocity curve, which we will obtain 
in the future. 

''Most of the uncertainty for F850LP is due to the aperture correction, so it should not be applied to 
magnitude differences within the same filter. 
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For each of the stars we have fixed the va lues of the teraperatu re and gravity derived from 
their spectral types and the calibration of iMartins et al.l (120051 ) and we have left -R5495 and 
£'(4405—5495), the monoch romatic equiva l ents o f Ry and E[B — V), as free parameters using 
the extinction law family of ICardelli et al. Jl989h. For each star we have used as input SEDs 
both t he CMFGEN grid calcu lated bv lMartins et all J2005 h and the TLUSTY OSTAR2002 
grid of iLanz fc Hubenyl (120031 ): however, after running CHORIZOS with the two SED fam- 
ilies, we could find no significant differences between them (not surprising, considering that 
they both have very similar intrinsic colors for a given temperature and gravity). Therefore, 
here we report only on the results for the CMFGEN case. Both Pismis 24-1 and Pismis 24- 
17 have 2MASS JHKg photometry, though Pismis 24-1 is unresolved. Given the increased 
accuracy of the extinction corrections when NIR photometry is included in CHORIZOS, we 
have done two different runs for each object: one with F330W-I-F435W+F550M+F850LP 
and another one with those same filters plus JHKg. Given that Pismis 24-1 is unresolved in 
2MASS, for its run with 2MASS we merged the HRC photometry for 24-lNE and 24-lSW. 

The results of the five CHORIZOS runs are shown in Table [2) the total extinction in 



F550M, the extinction law, the absolute visual magnitude (adopting from lMassey et al.ll2001 
a distance modulus of (m — M)o = 12.03±0.14), the apparent bolometric magnitude (i.e. the 
observed F550M magnitude with its extinction and bolometric corrections), the minimum 
X^, and the degrees of freedom of the fit. In all cases, a single well-defined solution is present 
in the likelihood i?5495 vs. £(4405 — 5495) plots. The mode SED is shown in Fig. [3] for 
Pismis 24-lNE (without 2MASS) and for Pismis 24-1 (with 2MASS). 



All the runs yield extinctions in F550M close to 6 magnitudes an d .^5435 values close t o 
3.0. The latter are significantly lower than the results found by iBohigas et al.l (120041 ). 



• Pismis 24-17 is m ore extinguished than Pismis 24-1 by ^0.4 magnitudes, which is 
consistent with the iBohigas et al.l (120041 ) results, though our values are lower overall. 



Table 2. CHORIZOS results. 



Object 


2MASS 


^F550M 


-R5495 


Mv 


{m- A + -BC)f550M 


Amin 


dof 


Pismis 24-lNE 


no 


5.54±0.10 


2.87±0.05 


-6.41±0.17 


1.71±0.10 


9.1 


1 


Pismis 24-lSW 


no 


5.52±0.10 


2.87±0.05 


-6.28±0.17 


1.73±0.10 


7.0 


1 


Pismis 24-17 


no 


5.91±0.11 


2.94±0.05 


-6.01±0.17 


1.91±0.10 


4.6 


1 


Pismis 24-1 


yes 


5.87±0.02 


3.01±0.01 


-7.50±0.14 


0.57±0.02 


40.4 


4 


Pismis 24-17 


yes 


6.34±0.02 


3.11±0.02 


-6.50±0.14 


1.42±0.02 


26.1 


4 
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3.0 2.5 2.0 1.5 1.0 0.5 



Fig. 3. — Mode SEDs and synthetic photometry for two of the CHORIZOS runs. The 
observed photometry is also shown, with the horizontal error bars marking the approximate 
extent of each of the seven filters (from left to right, HRC F330W, F435W, F550M, F658N, 
and F850LP; 2MASS J, H, and Kg) and the vertical error bars (small in all cases) the 
photometric uncertainties. The F658N (Ha) filter is not included in the CHORIZOS runs 
but is shown here for comparison: note how the emission excess is larger for Pismis 24-lNE 
than for Pismis 24-1. 
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• The inclusion of the 2MASS data introduces sginificant changes in the results, increas- 
ing the measured F550M extinction by ^0.4 magnitudes and -R5495 by R:i0.15. 



The minimum ^ values are too high. We attribute this to the lCardelli et al.l (119891) ex- 
tinct ion law, which was derived with an older photometric calibration (see lMaiz Apellaniz 



20061 ) and using only moderately-extinguished stars. It also uses a rather unphysical 



seventh-degree polynomial in the optical region, which is responsible for the long- 
period oscillations in wavelength visible in the reddened SEDs in Fig. [3], an effect that 
is unlikely to be real and that artificially increases the values. 



We have used the CHORIZOS results to derive new evolutionary zero-a ge masses for the 
three objects in our sample. The evolutionary tracks have been obtained from lLejeune fc Schaerer 
( I2OOII ). Additional uncertainties of 0.1 magnitudes have been added in quadrature to the 
bolometric magnitudes to take into account the problems with the extinction correction and 
the (small) uncertainty in the effective temperatures. The 2MASS information was included 
in the calculations because its presence significantly reduces the systematic errors introduced 
by the seventh-degree polynomial used for the extinction law in the optical (in any case, it 
is always good practice to add NIR data to optical photometry when analyzing objects with 
moderate or high extinction). We obtain nearly identical values for the three objects: 97±10 
M© (Pismis 24-lNE), 96±10 M© (Pismis 24-lSW), and 92±9 (Pismis 24-17). The real 
mass for the unresolved components in Pismis 24-lNE must be, of course, smaller: if the 
luminosity is equally split and the temperature is the same for both, the resulting values are 
of the order of 64 M© for each. 

As a test, we can also compute the mass of Pismis 24-1 as if it were a single star (i.e. 
the W02 hypothesis) and we obtaiE§| 155±19 Mq. This value is much lower than those 
obtained by W02 (210 or 291 M©). Given that our value for My is very similar to th eirs, the 
difference must arise from the d ifferent temperatu re scales used ( iMartins et al.ll2005l . derived 
from line-blanketed models, vs. IVacca et al.lll996l ). The reduction in temperature for a given 
spectral type motivated by the change moves the star in an HR diagram towards the lower 
right (the vertical effect is due to the smaller bolometric correction). Such a displacement 
is almost perpendicular to the evolutionary tracks for O stars that are just evolving off the 
main sequence and, thus, is a very effective way of reducing the mass estimate. 

What else could be affecting the masses? Even though Pismis 24 suffers moderate ex- 
tinction {^Ay 6 magnitudes), future improvements in this respect are likely not to be too 



^Note that the evolutionary track with the highest mass is 120 Mq, so we are forced to extrapolate, a 
procedure that should introduce additional uncertainties. 
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large, since our analysis includes multiband optical and NIR data simultaneously and Ak 
is only ^0.6 magnitudes. Also, our HST optical photometry is compatible with modern 
ground-based results and Tycho did not detect photometric variability, so the possibility 
of undetected changes (other th an those caused by t he known eclipses) seems unlikely. Fi- 
nally, the distance calculated by iMassey et al.l (120011 ) could be wrong: we plan to test that 
possibility in the near future by making use of our additional HST data. 



4. Conclusions 



We have presented HST and ground-based data that clearly resolve Pismis 24-1 into 
two objects separated by 363.86 mas and that indicate the existence of at least a third 
unresolved component in Pismis 24-lNE. We have also derived zero-age evolutionary masses 
for the three most massive objects at the core of Pismis 24-1, with values of 92 to 97 M©, 
one of them corresponding to the unresolved system Pi smis 24-lNE. The se values are very 
large, making the cores of Pismis 24-1 and Trumpler 14 (INelan et al.l 120041 ) the two locations 
within 3 kpc of the Sun with the highest density of very massive stars. However, these 
masses are not above what it is currently thought to be the stellar upper mass limit. 



After this manuscript was submitted, Virpi Niemela, one of the authors, passed away. 
The rest of the authors want to dedicate this article to her memory. She was the first one of us 
to study Pismis 24-1 and also the first one to suspect its multiplicity. But, more importantly, 
she was a wonderful human being to whom a generation of astronomers is indebted for her 
insight and her warmth. 

We would like to thank Phil Massey and his team of collaborators for letting us know 
about the optical photometric variability prior to publication and David Osip and Mark 
Phillips for kindly allowing us to use IMACS during part of an engineering night at the 
Baade telescope. We also acknowledge the useful comments made by an anonymous referee. 
Support for this work was provided by (a) the Spanish Government through grant AYA2004- 
08260-C03; (b) NASA through grant GO-10602 from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy Inc., under 
NASA contract NAS 5-26555; and (c) lALP-CONICET, Argentina. 
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